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Abstract 

Though topological aspects of energy bands are known to play a key role in quantum transport 
in solid-state systems, the implications of Floquet band topology for transport in momentum space 
(i.e., acceleration) are not explored so far. Using a ratchet accelerator model inspired by existing 
cold-atom experiments, here we characterize a class of extended Floquet bands of one-dimensional 
driven quantum systems by Chern numbers, reveal topological phase transitions therein, and the- 
oretically predict the quantization of adiabatic transport in momentum space. Numerical results 
confirm our theory and indicate the feasibility of experimental studies. 

PACS numbers: 05.60.Gg, 03.75.-b, 05.45.-a, 03.65.Vf 
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In both classical mechanics and quantum mechanics, position and momentum variables 
form a conjugate pair and can hence be treated on the same footing from a phase space 
perspective. The real physical world, however, does not have position-momentum symmetry. 
For example, energy bands of a solid are formed because electronic Hamiltonians are periodic 
in position but not in momentum. Due to such unequal roles of position and momentum, 
the mapping of quantum transport phenomena from position space to momentum space 
is typically nontrivial but, where possible, may lead to important insights and unforeseen 
opportunities. Anderson localization, for instance, was first discovered as a seminal result 
of quantum transport in position space. Its analog in momentum space was later found 
to be behind the intriguing phenomenon of "dynamical localization" [TJ [2] . This mapping 
has stimulated fruitful studies of Anderson transition in driven cold-atom systems [3j. As 
a second example, ratchet transport, namely, directed transport under a zero mean force 
in position space, has been mapped to momentum space as well, leading to the finding of 
ratchet accelerators (RA) IIHO]. 

Energy-band topology is of fundamental interest to studies of quantum transport in 
condensed-matter physics. The issue to be addressed here is whether transport in momentum 
space (i.e., acceleration) can be connected with band structures due to momentum space 
periodicity. At first glance this sounds impossible because, with the (non-relativistic) kinetic 
energy being a quadratic function of momentum, a realistic Hamiltonian is never a periodic 
function of momentum. However, for systems periodically driven by impulsive fields, the 
Floquet operator can still be periodic in momentum, thus forming Floquet (quasi-energy) 
bands. The topology of such Floquet bands then provides a useful tool in characterizing 
topological phase transitions in driven quantum systems [7] . How Floquet-band topology is 
manifested in acceleration then becomes an intriguing question. 

Using a one- dimensional RA model inspired by existing cold-atom experiments, we show 
in this work that the Floquet bands, defined on a 2-torus (formed by one Bloch phase and 
one experimentally tunable parameter) may be characterized by Chern numbers. We then 
theoretically show that adiabatic transport in momentum space can be quantized according 
to these topological numbers. Unlike quantized adiabatic pumping in position space [8], 
there does not exist a general (system-independent) flux operator in momentum space and 
hence the found quantization is about a quantized net change of momentum expectation 
value, rather than a pumping of particles through a cross section. Finally, though adia- 
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batic manipulation of Floquet states is often considered to be subtle [9j HP] , our numerical 
results confirm our theory and indicate that quantized momentum space transport can be 
observed in a wide regime, by simply scanning one system parameter in a relatively small 
number of discretized steps. Breakdown of quantization in momentum space transport may 
then be considered as a diagnostic tool for detecting non-adiabatic effects in Floquet state 
manipulation. 

To physically realize the main physics we need three ingredients: momentum space peri- 
odicity, an experimentally tunable periodic parameter and well-gapped Floquet bands. To 
be specific we consider a RA model [TT], which is based on an atom-optics realization of 
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a double- kicked rotor system [12]. The RA Hamiltonian is given by H = y + K cos(g + 
a) J2 n 8(t — nT) + K cos(g) J2 n ^(t~ nT — T ), with a corresponding Floquet propagator [11] 
U(a) = e -i(T-T )( P y2h) e -i(K/h)cos( q ) e -iT (p 2 /2h) e -i(K/h)cos( q +a)^ where all quant ities are properly 

scaled and hence in dimensionless units, q and p are canonical coordinate and momentum 
operators for cold atoms, and the 5 kicks are of period T, experimentally implemented by 
two optical-lattice potentials mutually phase-shifted by a, with equal strength K, and a 
time lag Tq <T [12]. Note that later a will be adiabatically tuned. Because the potential 
function is periodic in q, momentum eigenstates take eigenvalues mh + (3h, where G [0, 1] 
is a conserved quasi-momentum variable and m is an integer. To yield a Floquet opera- 
tor periodic in momentum despite the p 2 /2 term in the Hamiltonian, we set /3 = 0, which 
may be approximately implemented by considering a Bose-Einstein condensate whose co- 
herence width spans across many optical lattice constants [6j [13j E] . Effects of nonzero (3 
values will be discussed in Appendix [15J. If we now impose the quantum resonance con- 
dition Th = Att that has been one experimental subject [6j [13j HH IT6HT8] . we obtain an 
"on- resonance double-kicked-rotor model" [19J, with the Floquet propagator reducing to 



U R (a) = e^tf e-^^^e^tf e-^ cos(9+a) , (1) 

where h e = hT is an effective Planck constant, K e = KT , and a rescaled momentum oper- 
ator p e = T p. From here on, momentum exclusively refers to p e and we denote momentum 
eigenstates by |m), which has eigenvalue mh e and is periodic in q with period 2ir. Equation 
(10) indicates that if h e = 2itM/N, with M and iV being integers, then the Floquet operator 
U R is perfectly periodic in momentum space with a period of Nh e . According to Bloch's 
theorem, this momentum space periodicity leads to Floquet bands. Indeed, in the case of 



a = 0, the corresponding Floquet band structure [19] resembles Hofstadter's butterfly [20J, 
with h e here identified as an analog of the magnetic flux. It is such a remarkable resemblance 
(which also hints desirable band-gap features) that motivated us to connect the topological 
aspects of the Floquet bands with momentum-space transport. 




FIG. 1. (color online) Floquet eigenphase uj n (4>, a) vs and a for (a) K e = 3h e and (b) K e = 4h e . 

The Floquet band structure may be characterized by topological Chern numbers, provided 
that the bands are defined on a 2-torus. To that end, we supplement the Bloch phase 
parameter in momentum space with the periodic parameter a. This procedure somewhat 
resonates with recent efforts in identifying analogs of quantum Hall effect in one-dimensional 
systems [21]. The eigenstate-eigenvalue problem for U R now becomes Ur(o) \i[) n ((f),a)) = 
e iw„(ct>,a) |-0 n (0 ) q;^ 5 where G [0,27r) is the Bloch phase in momentum space and u n ((f), a) 
is the eigenphase of Ur{o). For a fixed pair of a and 0, N eigenphases can be obtained, 
and scanning (0, a) over [0, 27r) x [0, 2ir) forms extended Floquet bands on 2-torus, so-named 
to distinguish them from the common bands involving only the Bloch phase parameter 0. 
As a result n (1 < n < N) becomes an index of such bands (with h e = 2ttM/N). The 
eigenstates \ip n (<f>, a)) are chosen such that they are locally single- valued functions of a and 
0, and periodic functions of 0. Two computational examples of the Floquet bands o> n (0, a) 
are depicted in Fig. [T] for N = 3. It is seen that as the driving strength K varies, the 
landscape of the bands changes. Other calculations show that if h e = 2nM/N with even 
N, then there will be two (extended) bands touching each other. For simplicity, here we 
consider only odd N so that only accidental band collisions occur. 

In Fig. [7] we present the Chern numbers (see its surface integral expression below and 
see [15] for computational details) for a 3-band case. It is seen that at some isolated critical 
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FIG. 2. The Chern numbers of the 3 bands vs K e /h e , for h e = 2n/3. For results of a 7-band case, 
see Appendix [15] . 

values of K e , the Chern numbers jump, signaling the presence of toplogical phase transitions 
in our RA model. Indeed, the Chern numbers are invariant integers with respect to smooth 
deformation of the bands and only change discontinuously due to band collisions. Note that 
mainly in the context of quantum-classical correspondence in classically chaotic systems, 
Ref. |22j (see also Ref. [23]) formally studied the topological aspects of Floquet bands defined 
on a 2-torus formed by Bloch phases in both position space and momentum space. Our study 
is much different because (i) this work is based on an explicit physical implementation of 
momentum-space periodicity, (ii) here the Bloch phase in position space is fixed at (3 = in 
theory (so as to obtain Hofstadter's butterfly Floquet spectrum), and (iii) our Floquet bands 
are defined on a 2-torus that involves one experimental parameter, a key starting point for 
our theory below. 

We now seek the implications of the Floquet-band Chern numbers by considering an 
adiabatic cycle during which a increases from to 2tc. We first construct an initial state of 
the following form, 

\* n {a = 0)) = -^J rf0|^(0,a = O)), (2) 

which represents an equal-weight superposition of all the Floquet eigenstates of band n 
with a = 0. This coherent superposition state, which can be interpreted as a Wannier 
function in momentum space, uniformly samples all the Bloch eigenstates with different 
values of <fi (but all with a = 0), with a profile localized in the momentum space. Indeed, 
each eigenstate \if) n ((/), a. = 0)) is infinitely extended in momentum space, but |\l/ n (a! = 0)) is 
normalized to unity (localized in momentum space with Gaussian-like tails in all the cases 
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we studied). It is worth noting that because each eigenstate \ip n (4>,a = 0)) is defined only 
up to a global phase, one is free to choose an overall phase convention of \ip n ((j),a)) such 
that the superposition state in momentum space tends to be well-localized, thus making 
experimental preparation of the initial state easier. Such states can be highly localized so 
long as K is not too large. For example, the shown superposition state in Fig. 3(a) mainly 
occupies one momentum eigenstate, with small weights distributed over only a few nearby 
components. Given previous experiments where momentum superposition states in the same 
context were prepared [BJ, states as shown in Fig. 3(a) should be reachable in experiments. 
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FIG. 3. Momentum distribution \{m\^ 3 (a = 0))| = |^(m)| 2 [see Eq. (2)] with K e = 2h e for t = 
in (a) and after a 100-period adiabatic cycle in (b). Numbers shown in (a) are the phases of each 
momentum component. 



Consider then an adiabatic change in a, through a discretized protocol a s = 2-Ks/sf for 
(s — 1)T < t < sT, so that a completes one adiabatic cycle at t = SfT [24]. Assuming 
adiabatic following of the Floquet states, the state evolved from |^ n (a = 0)) [see Eq. (|2])] 
should stay as a superposition state at t = (sT)~ , with each component still being the 
eigenstate of Ur(o), with a = a s . That is, at t — (sT)~ , under adiabatic approximation the 
associated time-evolving state of the system becomes 

r>2lT 



|*n(a»)) 



1 

2^ 



dcj>\ilj n {^a s ))e ie ^\ 



(3) 



where 8 ((f), a s ) is the sum of a dynamical phase and a geometrical phase accumulated by the 
component starting from a = 0)). As a consequence of choosing \4> n ((fi, & = 0)) to be 

a periodic function of 0, 9(<p, a s ) is also necessarily periodic in 0. 

Next we evaluate (p e (s)) = (ty n (a s )\p e \ty n (a s )), namely, the momentum expectation 
value on the state \^f n (a s )). To proceed we first write Floquet eigenstates \ip n (<j>, a)) as the 



product of two parts using the Bloch theorem, namely, \^ n ((j),a)) = X((p) \u n (4>,a)), with 
X{4>) = e m ^ Nhe . We then define \u n ((f),a)) = J2m=i \ m )( m \u n {4>, a)). It can be shown 
that the iV-element state \u n (<p, a)) is an eigenstate of the following N x N reduced Floquet 
matrix 

oo 

(m\U(4>, a)\m') = (m|^W^(«)^^h' + IN), (4) 

l=— oo 

with m, m! — 1, • • • , N. Lengthy but straightforward calculations [J5] then yield a compact 
expression for (p e ), i.e., (p e ) = ^ J Q 27r dcj){u n ((f), a s )\ih e -^\u n ((j), a s )). To compare this expec- 
tation value with its preceding value for a = a s _i, we consider a first-order perturbation 
theory so as to rewrite \u n (<f), a s )) in terms of \u n (<p, a s _i)}. Specifically, to the first-order of 



5a = a s — a s _i = 2ix/sf, we have [15 



\u n ((j),a s )) = \u n {(f), a s _i)) 



g«w n (<j!),o s _i) _ giaj n /(</),a s _i) 
n'=l,y^n 

x |tt n ,(0,a s _i)). (5) 

The change in momentum expectation value (denoted by 5{p e ) s ) over the period from t = 
(s — 1)T to (sT)~ can now be calculated to the first order of 5a, 

1 f 2n r d 

5(p e ) s = —N J d(p^(u n ((/),a s )\ih e — \u n ((f),a s )) 

- (u n ((f), a s _i)|zfr e — \u n {(j), a s _i)) k (6) 



Further substituting Eq. (22) into Eq. (21) yields 

r-2?i 



N f 

$(Pe)s = -— d(f) B n ((f),a s ) 5a, (7) 
27T y 



where B n ((f), a) is identified as the Berry curvature 



B n ((j),a)=i ^ 



^ {u n \ aj. \u n ') (u n i | a \u n ) 



(8) 



with the explicit dependences of \u n {4>, a)), u n ((p, a) and U((p, a) on a and all suppressed 
for brevity. 



The total change in momentum expectation value from t = to (sT) is denoted by 
A Pe (s). Because A Pe (s) = Y? a '=i HPe)s', we have 

A Pe (s/) = -M/ d<f> dctB n ((f>,ct) 
Jo Jo 

= -2TrMC n , (9) 

where we have used H e = 2irM/N and C n is exactly the Chern number of the nth Flo- 
quet band, whose surface integral expression is C n = j- d<p da B n ((J),a) [15J. Thus, 



Eq. (28) reveals that the net change in the momentum expectation value over one adiabatic 
cycle of a (starting from state ^^(a = 0))) is quantized: it should be proportional to the 
Chern number of the nth extended Floquet band. This is our central theoretical result. 

It is necessary to numerically verify our theoretical insights above. Detailed results are 
shown in Fig. [ij again for the case of h e = 2n/3, with the adiabatic cycle lasting for Sf = 100 
periods (in this case, see also Fig. 3(b) for the final momentum space profile) or Sf = 1000 
periods. First of all, apart from the regime of K e values near the critical point K e pa 4.2h e 
(see Fig. 2), our numerical values of A Pe (s/)/(— 2n) almost perfectly match the Chern 
numbers. This is the case before or after the jumps of the Chern numbers. The insets of 
both panels depict how A Pe (s)/(— 2tc) builds up with time and eventually reaches integer 
values that match the Chern numbers. In addition, we have checked numerically that if 
we repeat the adiabatic cycle, then the same quantized increase in momentum expectation 
value is obtained (15] . We are thus witnessing a clear quantization effect in acceleration as an 
outcome of Floquet band topology. Note however, in the vicinity of phase transition points, 
e.g., K e /h e pa 4.2, momentum space transport is no longer quantized. This is because if a 
topological phase transition is about to occur, then the associated band gaps are not large 
enough to guarantee adiabaticity. Supporting this understanding, a comparison between 
Fig. [4] (b) and Fig. [41(a) shows that a longer adiabatic cycle indeed significantly narrows 
down the non-quantization window. Our numerical data suggests that at least for the 3- 
band case here, if the driving field strength K e is far away from the phase transition points, 
then only 50-100 kicking periods (depending on K e ) will be needed to observe quantized 
acceleration. This is experimentally motivating, because Floquet state manipulation itself 
is a topic of much theoretical interest (9J \10\ ■ The robustness of this quantization against 
perturbations is also examined in Appendix [IB] . In short, for the 3-band case above, the 
quantization effect can tolerate about 0.5% uncertainty in h e and a nonzero (3 around 0.01, 
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FIG. 4. Change in the momentum expectation value (divided by — 2ir) vs K e /h e , after one adiabatic 
cycle implemented in (a) 100 and (b) 1000 discretized steps, for initial states prepared on each of 
the three Floquet bands [see Eq. Q]. Insets shows A Pe (s)/(— 2tt) vs number of periods s, for 
K e /h e = 2.0 (K e /h e = 6.0) on the left (right). In each inset, each of the three plotted curves is for 
one of the three Floquet bands, which in the end approaches integer values that match the Chern 
numbers. 

which should be achievable in light of previous experiments [T3] . 

In conclusion, based on extended Floquet bands, we have exposed topological phase 
transitions in driven quantum systems and demonstrated how quantized adiabatic transport 
in momentum space may emerge from Floquet band topology. Numerical results based on 
a cold-atom-based dynamical model suggest that future experimental verification of our 
results is possible in terms of initial state preparation, adiabatic cycle implementation, and 
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the robustness of quantization. 
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APPENDIX A: DERIVATION OF QUANTIZED MOMENTUM TRANSPORT 



We use the same notation as introduced in the main text. The Floquet propagator of a 
on-resonance double-kicked rotor is given by 

with eigenstates and eigenvalues denoted by 

U R (a) \M<f>,<*)) = e^ (0,a) lV>n(<M) • (11) 

Here \ip n (<p,a)) represent the Bloch eigenstates of Ur(o) with Bloch phase 0, due to 
momentum-space periodicity of period Nh e for h e = 2nM/N. Using the Bloch theorem, 
one rewrites (m \ tp n ((fi, «)) as 

(m <*)> = e im<t> ^(m \u n (<f>, a)) . (12) 

Here (m + N\u n (<f), a)) = (m\u n (4>,a)) holds for all momentum eigenstates \m). We adopt 
the normalization convention Xlm=i I ( m \ u n(4>y a )) | 2 = 1- With this normalization conven- 
tion, it can be shown that the identity operator in terms of the Bloch eigenstates is given 
by 



-27T N 

27T 

and their inner products are given by 



- #V|^,a))(^,a)| = l, (13) 



(M<t>')\M<f>)) = 2vr5„,„^(0 - 0'), (14) 



which implies that the initial states we considered in the main text [see also Eq. ( 15 ) below] 
are normalized to unity. Note that we have chosen a basis of Bloch eigenstates, {\ip n (<p, a))}, 
which are set to be periodic functions of 0, meaning that \if) n {<f> — 0,a)) = |"0n(0 = 27r,a)) 
for all a. Also, the overall phase convention of {\ip n (<f),a))} can be chosen such that it is a 
locally single- valued function of and a. 

Next, we consider the states involved in the adiabatic cycle described in the main text. 
Such states are of the form 

2tt 



l*n(«.)> = ^ / ^0 l^n(0, a.)) e w ^\ (15) 



12 



where s takes values from to Sf as explained in the main text and 9((j),a = 0) = 0. 
We begin our derivation of quantized transport in momentum space by taking momentum 
expectation value of such states, 

oo 

(V n (a s )\p e \V n (a s )) = mh e (m\-$ n (a s )) (V n (a 8 )\m) . (16) 



m=— oo 



We insert two copies of Eq. (13), one in the middle of each of the inner products of the right 



hand side. Then, making use of Eqs. (15), (14) and (12), we arrive at 



El If 17TL\<p — (J> ) j 

mh e / d<p I dcfie * (m\u n ((j), a s )) (u n (<p , a s )\m) 

m=-oo J ° J( > 

x 1 ^ c i(9(J>,a)-e(J>',a))_ ^ 

We then rewrite m as a derivative in </> acting on the exponential term and perform an inte- 
gration by parts to shift this derivative off the exponential and onto the other (^dependent 
terms. Further, we split m into m' + sN, where m' takes integer values 1,2, ••• ,N, 
and s takes all integer values from — oo to oo. Then, after making use of the fact that 
(m + sN\u n ((p, a)) = (m\u n (<j), a)) for all momenta mh e and all s e Z, one invokes Poisson's 
summation formula Yl^L-oo e ls ^~^> = Yl^L-oo 2tt<5(0 — (j)' — s 2ix) and finds 

2~ 



1 N f 2n d 

2^ z\2 J d(j)(u n ((p,a s )\m') Nih e — ((m'\u n ((f),a s ))) 



m'=l ' 

2, Q 

d(f) (u n ((f), a s ) | Nih e — \u n ((f), a s )) , (18) 



o 



2tt 

where \u n (4>, a)) = Ylm=i \ m ) { m \ u n{.4>i a )) an d we have made use of the fact that dej) 99 ^?*' 
0. Summarizing, we have shown that 

1 f 2w d 

(^ n (a s )\p e \ty n (a s )) = — / d(j)(u n ((j),a s )\Nih e —\u n ((j>,a s )). (19) 



o 



2tt 

Note that \u n ((f>,a s ) are the eigenstates of the reduced Floquet matrix 



(m\U(<t>,a 3 )\m') = (m\Xi(<f>)U R (a a )X((i>)\7ri + IN), (20) 

l=— oo 

with m, m' = 1, • • ■ , N and X(4>) = e ip ^ /Nhe . 

Next, we consider the difference in momentum expectation value between consecutive 
states in the adiabatic cycle, namely, 

5( Pe ) s = (^ n (a s )\p e |* n (cO) - (* n (a a _ x )| J) e |^nK-l)) ■ (21) 
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For both terms in the above equation, the expectation values are given by Eq. (19) with 
a = a s _i and a = a s respectively. Next one may express \u n (<f>, a s )) in terms of \u n (<p, a s _i)} 
to the first order in 5a = a s — a s -%, i.e., 

\u n {(f), a,)) = \un(<f>, a s -i)) + \u£\(i), a s _i)) (22) 

where 

" (u n ,(<f), a a - 1 )\ dD ^' as - l) \un(<t>, a B _i)) 
\u^(<f>, = 6a ^ e ^ n (0,^-i) _" e W(^-i) M&a-i))- (23) 

n'=l,^n 

This perturbation theory can be performed because (i) our Bloch eigenstates can be chosen 
as continuous in a for any fixed <j) an d (ii) in the adiabatic limit Sf — > +oo, 5a — > and 
hence only first-order matters. It is now seen that \u n ((f), a s )) is expressed as the eigenstate 
(correct to the first order in 5a) of the reduced Floquet matrix with a = a s _i, plus a 
correction. Note also that the expression dU Q^ as ~^ in the correction term can be explicitly 
calculated by U a via 

(mlUfaa^m!) = V (m\X\<j>) f^l x{(j>)\m' + IN). (24) 

t = — OO 

Finally, we use the Feynman-Hellmann-like formula 

d \ (-u w (0, a)| ^-M^ \u n '(4>, a)) 
-(u n {^a)\ \ \u n ,(<P,a)) = — _ — , (25) 



and its adjoint. Substituting Eq. (25), its adjoint, and Eq. (22) into Eq. (21) and keeping 



only first order terms in 5a, we obtain 

5( Pe )s = -^~Nh e [ d4> B n (<f>, a s ) 5a, (26) 
2tt J 



where B n ((j),a) is identified as the Berry curvature 



n'=l,^n I 1 

with the explicit dependence of \u n ((f),a)), uj n (<p,a) and U(<f),a) on a and <p all suppressed 
for brevity. 

The total change in the expectation value of p e from t — to (sT)~ is denoted by 
A Pe (s). Because A Pe (s) = X/s'=i ${Pe)s', one finds the total momentum change over the 
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entire adiabatic cycle, i.e., 

p2tt p2tv 

A Pe(sf) = -M d(j) daB n ((f),a) 
Jo Jo 

= -2irMC n , (28) 

where we have used h e = 2nM/N and C n is exactly the Chern number of the nth Floquet 
band defined earlier, whose equivalent expression is 

/*2"7T /•27T 

Cn = 2^J d(l> J daB ^ a )- ( 29 ) 
This completes our derivation of quantization of transport in momentum space. 



APPENDIX B: DETAILS OF CHERN NUMBER CALCULATION 

Here we present some details regarding how the Chern numbers of Fig. 2 of the main 
text are calculated. In terms of a line integral, the Chern number is given by 

c n = ±f@ n (r)\-^\j n (r))dr, (30) 

where T = (0, a), |i/> n (r)) = Xlm=i \ m ) ( m lV'n(r)) (normalized on N momentum compo- 
nents). We first mention that this expression is equivalent to ^ § (u n {Y)\ \u n {T)) dT , 
where \u n (T)) = E™=i l m ) ( m M r ))> and |^n( r )> = X(<f)) \u n (<f),a)), with X(<f>) = 
e %p e (t>/Nhe_ rp see this, one needs only to substitute \ip n (T)) = X(<f>)\u n (</>,a)) into the 
Chern number expression above and write the integral as the sum of 4 terms, each on 
one edge of the (0, a) Brillouin zone. The two terms containing S- immediately reduce 
to the form ^ J (u n (T) \ \u n (T)) da because X((p) is not a function of a. The other two 
terms containing will also reduce to the similar form ^ J {u n (T)\ \u n (T)} d<p after one 
application of the product rule of differentiation and noting that two of the resulting four 
terms cancel. This then leaves us with C n = ^ § (u n (T) \ ^ \u n (T)) dT. We can now work 
on C n in terms of \u n ) to make the link with the main text more concrete, though it should 
be clear that it is entirely valid to work with \ijjn) instead. In the remainder of this section, 
we will derive the surface integral formula for the Chern number as it was used in the main 
text and then mention some details on our numerical calculation. 

The surface integral for C n which involves the Berry curvature as shown in the main text 
is obtained from the line integral via the steps shown here. First, a single application of 
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Stoke's theorem leaves us with 

r-2-K r2ir 



C n 



2ty 

% 

2^ 



d(f)da 



o jo 

2tt r2ir 



d ( , d \ d ( , d 

(u n ((f),a)\ — \u n ((f),a)) )-Q^[ (u n ((p,a)\ ^- \u n ((f),a)) 



o Jo 



d(f)da 



d \ d 

— (u n ((f),a)\ 1 — \u n ((f),a)) 



d \ d 

— (u n ((f),a)\J — \u n ((f>,a)) 



(31) 



We then insert one copy of the N x N identity matrix of the form J2 n ' (^n'l — ^-nxn 
between each of the two inner products, where we have refrained from writing the functional 
dependence of \u n ) on (0, a) for brevity. The term with n! = n in one inner product 
will cancel that in the other inner product because (j^ (u n \ \ \u n ) = — (u n \ (j^ \u n )^j and 
similarly for This leads to 



2tt 



2n r 2n 



JO 



N 

d(f>da 



l~{u n \) \u n ,)(u n ,\^\u n ) 



^ (u n ((f),a)\^ \u n >) (u n ,\ ^ 

(32) 

To proceed, we note that the states \u n (<p, a)) are eigenstates of the reduced Floquet matrix 

oo 

i\U{4>, a)\m') = (m\X\4>)U R (a)X{4>)\m' + IN), (33) 



(m 



l=— oo 



with m,m' = 1, • • • ,N, and where Ur(qz) is the "on-resonance double kicked rotor model" 
Floquet operator defined in the main text. From the eigenvalue equation U(<j), a) \u n ((f), a)) = 
e iui n ((f),a) \u n ((J), a)), we may take a derivative with respect to <fi on both sides, take the inner 
product with (u n '^ n ((f), a) \ from the left on both sides, and rearrange to obtain the Feynman- 
Hellmann-like formula 

v(0,a)| 9U ^ a) \u n {(f),a)) 



(u n/ (4>,a) \ — \u n ((f),a)) -- 



(34) 



which has adjoint given by 



— {u n [(j),a)\ \u n i{(f>,a)) = 



)-!W„ p — lU! I 



(35) 



Similar expressions also hold for the S- derivative. We substitute the Feynman-Hellman-like 



formulas and their adjoints for both and a into Eq. (31) and obtain 

»2tt r2-K 



C n 



2tt 



o Jo 



d(f)da B n ((f>, a), 



(36) 
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where 

n'=l,^n I 1 1 J 

is the Berry curvature. This completes the derivation of the surface integral for C n featured 
in the main text. 

We now mention some brief details on how we numerically calculate the Chern numbers. 



We chose Eq. (36) to calculate the Chern numbers rather than Eq. (30) because the former 
method allows us to use any phase convention for the \u n ) eigenstates (a quick check shows 
that the Berry curvature is invariant with respect to replacing each eigenstate with itself 
multiplied by a global phase). First, we discretize the (<ft, a) Brillouin zone into a sufficiently 



fine grid. We then numerically work out the integrand of Eq. (36) at each point on the grid 
and numerically integrate to obtain the Chern number. The different parts of the integrand 
are obtained as follows. We numerically evaluate the N x N elements of the U(<f>, a) matrix 



via Eq. (33). The infinite sum on the right hand side there does not present a problem 
because (m\Uii(a)\n) decays rapidly as \m — n\ increases, so we can truncate the sum at 
some suitable point with practically no loss of accuracy. Each term (m\Uii(a)\n) is efficiently 
calculated via Fast Fourier Transform (FFT). We then numerically diagonalize and obtain 
the iV-element vectors representing the \u n ) eigenstates as well as the eigenphases u n . Lastly, 
we obtain the elements of the N x N matrices |^ and This is done by acting -j^r and 



on both sides of Eq. (33), simplifying the expression on the right hand side by hand, and 
then obtaining the numerical values via FFT. The infinite sums here may again be truncated 
at a suitable level with no appreciable loss in accuracy. 



APPENDIX C: STABILITY OF QUANTIZED TRANSPORT TO PERTURBA- 
TIONS 

To motivate potential experimental interest, in this section we consider the effects of of 
nonzero quasimomentum variable /3, small perturbations to effective Planck constant h e , 
and inevitable imperfections in actual experimental implementations on the quantization of 
transport in momentum space. 

To consider the effects of nonzero /3, we begin from the double-kicked rotor model Floquet 

17 



propagator, 

(j( a ) = e - i ( T - T o)(pV^) e -KK/h)cos(g) e -iT (p 2 /2h) e -i(K/h)cos(q+a)^ /gg-j 

which describes a kicked particle moving on a line. As explained in [25], the periodicity of 
the Floquet propagator in q allows us to map the dynamics of a kicked particle on a line 
onto the dynamics of a fictitious ensemble of kicked particles each on a circle, with every 
such particle representing the dynamics for a different /3-component of the kicked particle 
on the line. Each such particle on the circle is thus referred to as a /3-rotor. The rotors 
may then be evolved separately from one another and recombined appropriately at the end 
to recover the time-evolved state of the actual particle on the line. The different /3-rotors 
evolve independently of one another under the action of a Floquet propagator 

Ub(oc) = e - i ( T - T o)((n+/3rOV2^) e -i(^/^)cos(5) e -iT ((n+/3fi)V2rO e -i(^/fi)cos(g+ a )^ /gg) 

where n + 0h is the /3-rotor's momentum operator with eigenstates \m) and momentum 
eigenvalues (n + 0)h with n G Z. Now, imposing the main quantum resonance condition 
TH = 47r, and writing h e = hT , K e = KT and h e = nT , the Floquet propagator becomes 

Up(a) = e - i ^- 1 )i^e+me) 2 /2H e ] e -i(K e /h e )co S {q) e -i[(n e +^ ( 4 q) 

The rescaled momentum eigenvalues are then given by (n + 0)h e , neZ. If we set = 0, the 
momentum operator has eigenvalues of only integer multiples of h e and we recover the on- 
resonance double-kicked-rotor-model Floquet propagator described in Eq. (1) of the main 



text. For nonzero (3, Up in Eq. (40) lacks translational invariance in momentum, so the 
entire derivation in the previous section will not apply and we hence do not expect to see 
quantized momentum transport in the /3-rotor dynamics. However, for small values, it 



is reasonable to expect that evolution of the special initial state in Eq. (15) will still show 
quantized transport up to a good approximation. We find that this is indeed the case, as 
evidenced in Fig. 5(a), where we see that for values up to 0.01, the quantized transport still 
survives for K e = 2.0h e , h e = 2n/3. Experimentally, however, the initial states prepared are 
typically a mixture of states with values following a narrow Gaussian distribution peaked 
at = 0. Hence, to gauge the experimental feasibility of our proposal, we computed the 
transport values for a range of values near and took the average of these values weighted 
with a Gaussian peaked at 0. We found that for a Gaussian distribution with standard 
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deviation of 0.01, the change in momentum expectation value divided by — 2ty for the 3 
bands at K e = 2h e for a 100-period adiabatic cycle are —0.99, 1.98 and —0.99 for bands 1,2 
and 3, respectively. In other words, the quantization of transport is still observable even 
with a realistic spread over (3 values. This indicates that the small {3 spread in current cold- 
atom experiments based on Bose-Einstein condensates should not hinder the observation of 
quantized transport in momentum space. 

Next, we briefly consider the effects of small deviations in the Planck constant H e for 
the same 3-band case considered above with h e = 2n/3 and K e = 2h e . Certainly, even a 
slight change in h e should cause the Floquet-band structure to change entirely (recall that 
the Floquet band structure depends on whether H e is rational or irrational), but physical 
observables should not be as sensitive as the band structure itself. Indeed, as seen in Fig. 
5(b), the quantized transport is still observable for deviations in h e up to 0.005ft e . 




FIG. 5. (a) Adiabatic momentum transport as a function of time for j3 = 0, 0.01, 0.02. (b) Adiabatic 
momentum transport as a function of time for 5H e = 0, 0.005/i e , 0.017i e . In both cases, we consider 
adiabatic cycles lasting for 100 kicking periods, with the initial state |*3(a = 0)) prepared on band 
3, for K e = 2.0h e and H e = 2ir/3. Quantized adiabatic transport is observed in those cases where 
the final values of A Pe are close to the integer value —1. 

Finally, we consider the effects of practical imperfections in experiments. To take into 
account the difficulty of superposing a large number of different momentum eigenstates in 
actual experiments, we consider a perturbed version of the initial state ^(ck = 0)} by set- 
ting to zero all probability amplitudes outside the range m G [—5, 5] and multiplying the 
resulting state by an appropriate constant for normalization. This does not alter the original 
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FIG. 6. Momentum transport as a function of time in the presence of noise (see the text for the 
noise details and system parameters). The case of panel (a) includes noise present in the amplitudes 
and phases of each momentum eigenstate component in the initial state. The case of panel (b) 
simulates noise present in the value of phase shift a during each step of one adiabatic cycle. In 
each case, we averaged over 1000 realizations of the noise. A and B represent noise intensity. The 
plotted error bars represent the standard deviation in the total momentum transport found in our 
numerical experiment. Note that the averaged total momentum transport stays very close to the 
quantized value despite the relatively strong noise. 

state significantly because its probability is tightly concentrated about m = (more than 
99.9 % of the original state's probability lies within m 6 [—5,5]) . With this truncated 
state as our starting point, we consider two noise models. In the first model, we introduce 
an uncertainty in the relative phases of the 11 surviving momentum states in the superpo- 
sition by multiplying each one by exp[z(0.05)27r£] where £ is a random variable uniformly 
distributed in [0,1). We then assume that the probability amplitude of each constituent 
momentum state is scaled by the term [1 + A(£ — 0.5)] and then multiply all states by an 
appropriate constant for normalization. In our second noise model, we assume that in the 
s-th period, the value of the phase shift a is given by [s + B(£ — 0.5)]27r/100. We now 
present in Fig. 6 two simple numerical examples for band 1 and band 3, with Sf = 100, 
K e /h e = 2.0, and h e = 2n/3. It is seen that in both models, the total momentum transport 
over one cycle (divided by — 2n), when averaged over the 1000 realizations of the random 
noise, changes only slightly from the quantized value (—1) despite the rather strong noise. 
The fair robustness to imperfection to initial state shown in Fig. 6(a) is partially due to the 
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fact that our initial state is dominated by very few momentum components. Indeed, the 
fluctuations in the total momentum transport from realization to realization, as manifested 
by the plotted error bars for band 3, are much smaller than those for band 1, consistent 
with the fact that our initial state for band 3 turns out to be more localized than that for 
band 1. Results in Fig. 6(b) illustrate that slow fluctuations in a (noise amplitude rather 
comparable to 2rr/sf) during the adiabatic process do not severely affect the quantization. 
As seen from Fig. 6(b), the total momentum transport averaged over 1000 noise histories 
is only slightly shifted from the quantized value (—1), and fluctuations represented by the 
error bars are also rather small. This confirms that the fine details of an adiabatic process 
are largely irrelevant. 



APPENDIX D: THE SEVEN BAND CASE 
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FIG. 7. The transitions of Chern numbers as kicking strength K e increases, for h e = 2tt/7. The 
horizontal lines represent the bands and the numbers written on them are their Chern numbers. 
Vertical lines connecting two horizontal lines represent degeneracies occurring between the two 
associated bands. Note the unusually big changes of Chern numbers across some critical points. 

In this section, we wish to emphasize that our topological characterization of the Floquet 
bands using Chern numbers is equally applicable to cases with more bands. Here we show 
a case with 7 bands, with h e = 2ir/7. We present in Fig. 7 the Chern number results as 
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the kicking strength increases. The associated Chern numbers are found to be symmetric 
about the central band, so we need only display 4 bands. The same qualitative features as 
we observed in the 3-band case are seen here. In particular, all bands other than the central 
band have Chern number —1 for low kicking strength, with the central band having a Chern 
number equal to the negative of the sum of the Chern numbers of all other bands. At band 
collisions, the Chern number of each band jumps by a multiple of N = 7 for h e = 2n/N. For 
this reason we see that the Chern numbers can rapidly become very large integers, which 
is a quite interesting feature for our system. We have also carried out adiabatic transport 
studies for this case. Due to the smaller band gaps for the 7-band case, we need to have 
longer a-cycles (i.e., a must be varied from to 2tt over a larger number of periods) in order 
for effective adiabatic following to take place. Experimentally speaking, this case is hence 
of less interest. However, in our numerical studies, other than the longer adiabatic cycles, 
the results are qualitatively the same as the 3-band case presented in the main text. 



APPENDIX E: FURTHER NUMERICAL DATA 



In this last section, we provide some further examples of distributions in momentum space 
after the adiabatic cycles proposed in the main text. We choose K e = 2H e for H e = 27r/3, as 
we did in Fig. 3 of the main text. Figure 8 shows the initial and final distributions after a 
single adiabatic cycle for initial Wannier states prepared on bands 1 and 2. Figure 9 shows 
the distribution in momentum space after completion of a second and third adiabatic cycle 
for an initial Wannier state prepared on band 3 (see Fig. 3 of the main text). 
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FIG. 8. For the initial state \^i(a = 0)) for K e = 2h e prepared on band 1, the distribution in 
momentum space before and after a 100-period adiabatic cycle are shown in (a) and (b) respectively. 
Panels (c) and (d) show the parallel results starting from the initial state | l I / 2( a = 0)) for K e = 2h e 
prepared on band 2. In panels (a) and (c), the numbers displayed are the phases of the amplitudes 
of the constituent momentum eigenstates in the initial superposition state. 
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FIG. 9. Final momentum distribution |(m|*)| after (a) two 100-period adiabatic cycles and (b) 
after three 100-period adiabatic cycles. The initial state was the Wannier state 1*3 (a = 0)} for 
K e = 2h e (same as that used in Fig. 3(a) of the main text). 
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